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Abstract. In this paper we present two independent computational proofs that the 
monoid derived from 5x5x3 contingency tables is normal, completing the classifi- 
cation by Hibi and Ohsugi. We show that Vlach's vector disproving normality for the 
monoid derived from 6x4x3 contingency tables is the unique minimal such vector 
up to symmetry. Finally, we compute the full Hilbert basis of the cone associated 
with the non-normal monoid of the semi-graphoid for \N\ = 5. The computations are 
based on extensions of the packages LattE-4ti2 and Normaliz. 



1. Introduction 

Let 5* = monoid (G) be an affme monoid generated by a finite set G C Z n of integer 
vectors. We call S normal if S = cone(G) fl lattice(G), where cone(G) = {x G W 1 : 
x = -^iSi) \ e ^+) Si e G} denotes the rational polyhedral cone generated by G 
and where lattice(67) = {x e I" : x = ^ Ajgj, Aj G Z, g { G G} denotes the sublattice 
of Z n generated by G. In this paper, we will stick to the case that lattice(G) = Z n . 
Then, normality of S is equivalent to saying that G contains the Hilbert basis of 
cone(G), i.e., every lattice point in cone(G) can be written as a nonnegative integer 
linear combination of elements in 67. By the Hilbert basis H{C) of a pointed rational 
cone C we mean the unique minimal system of generators of the monoid M of lattice 
points in C. The Hilbert basis of C consists of the irreducible elements of M, i.e., those 
elements of M that do not have a nontrivial representation as a sum of two elements 
of M (see [21 Ch. 2] for a comprehensive discussion). Note that deciding normality of 
an affme monoid is NP-hard [S]. 

Normality of monoids derived from r\ x r 2 x • • • x r^ contingency tables by taking 
N — 1-marginals (that is, line sums) was settled almost completely by Hibi and Ohsugi 
[H]. In this paper we close the last open cases by showing computationally, via two 
different approaches and independent implementations, that 5x5x3 has a normal 
monoid. The normality for 5x5x3 implies normality for the other two open cases 
5 x 4 x 3 and 4 x 4 x 3 by [9, 3.2]. 
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Here is the defining matrix A 5x5x3 whose columns generate the monoid associated to 
5x5x3 contingency tables. Every • corresponds to an entry 0. 
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Note that this normality problem cannot be settled directly by computing the Hilbert 
basis of the associated cone using state-of-the-art software such as Normaliz v2.2 [2111] 
or 4ti2 vl.3.2 [TJ [6]. Both codes fail to return an answer due to time and to mem- 
ory requirements of intermediate computations. Using the computational approaches 
presented below, we can now show the following. 

Lemma 1. The monoid derived from of 5 x 5 x 3 contingency tables by taking line 
sums (= two-marginals) is normal. 

This completes the normality classification of the monoids derived from n x r 2 x • • • x r^ 
contingency tables by taking line sums as given in [9]: 

Theorem 2. Let r\ > r2 > ■ ■ ■ > r^ > 2 be integer numbers. Then the monoid derived 
from ?"i x r 2 x • • • x r^ contingency tables by taking line sums is normal if and only if 
the contingency table is of size 

• r\ x r 2 , r\ x r 2 x 2 x . . . x 2, or 

• r\ x 3 x 3, or 

• 4x4x3, 5x4x3, or 5x5x3. 

For the monoid of 6 x 4 x 3 contingency tables, a vector disproving normality was 
presented by Vlach [12]. The right-hand side vector f for the counts along the coordinate 
axes is given by the following three matrices: 
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The unique point in the 6x4x3 transportation polytope {z e M 72 : Az = f , z > 0} is 
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So z* is indeed a hole of the 6x4x3 monoid. We are able to show the following. 

Lemma 3. The right-hand side vector f presented by Vlach [12] is the unique vector 
(up to the underlying SqX S^x S3 symmetry) in the Hilbert basis of the cone of 6x4x3 
contingency tables that is not an extreme ray. 

The treatment in [7] now completely describes all holes of the cone, that is, all lattice 
points in cone(A6 X 4x3) that cannot be written as a nonnegative linear integer combi- 
nation of the (integer) generators of the cone: 
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Corollary 4. Let f be the hole in cone(A 6X 4 X3 ) presented by Vlach [TjJj and let z* G M. 7 ^ 
be the unique solution to v4 6X 4 X3 z = f , z G , as stated above. Moreover, let G denote 
the set of those 24 columns of Aq X 4 X ^ for which z* > 0. 

Then the set of holes in cone(A 6x4x3 ) is the set of all points that can be written uniquely 
as er(f + s) with o G 5*6 X S4 x S3 and with s G monoid(G). 



Finally, we have computed the Hilbert basis of the cone associated to the semi-graphoid 
for \N\ = 5 [10J. It was already shown in [8] that the corresponding monoid is not 
normal by constructing a hole via a different method. The computation of the full 
Hilbert basis was not possible at that time, neither with Normaliz, nor with 4ti2. 
Here is the defining matrix whose columns generate the monoid associated to the semi- 
graphoid for \N\ = 5. Every . corresponds to an entry 0. + and — represent entries 1 
and —1. 
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Lemma 5. The Hilbert basis of the cone associated to the semi-graphoid for \ N\ = 5 has 
1300 elements that come into 21 orbits under the underlying symmetry group S$ x S2. 
These are represented by the 21 rows of the following matrix: 
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2. Computational Approaches 

In this section we present the two computational approaches that allowed us to solve 
the three challenging Hilbert basis computations of the cones associated to 5 x 5 x 3- 
tables, to 6 x 4 x 3-tables, and to semi-graphoids for \N\ = 5. In the first approach, we 
iteratively decompose the cone into smaller cones and exploit the underlying symmetry 
and set inclusion to avoid a lot of unnecessary computations. An implementation of 
this approach is freely available in the new release latte-f or-tea-too-1 .4 of "LattE 
for tea, too" (http://www.latte-4ti2.dej, a joint source code distribution of the two 
software packages LattE macchiato and 4ti2. In the second approach, we exploit the 
fact that the cones are nearly compressed; hence many cones in any pulling triangula- 
tion are unimodular, and the same holds in placing triangulations. Using our second 
approach, none of these unimodular cones is constructed, saving a lot of computation 
time. An implementation of this approach will be freely available in the next release of 
Normaliz ( ]http : / /www. math. uos . d e/norm aliz[ ), together with the input files of the 
examples of this paper. 

2.1. First approach: exploiting symmetry. Let us assume that we wish to compute 
the Hilbert basis of a rational polyhedral cone C = cone(ri, . . . , r s ) C R™. Moreover, 
assume that C has a coordinate-permuting symmetry group S, that is, if v e V and 
o G S then also er(v) e C. Herein, the vector <r(v) is obtained by permuting the 
components of v according to the permutation a . 
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One approach to find the Hilbert basis of C is to find a regular triangulation of C into 
simplicial cones Ci,...,C^ and to compute the Hilbert bases of the simplicial cones 
Ci, . . . , Ck- Clearly, the union of these Hilbert bases is a (typically non-minimal) system 
of generators of the monoid of lattice points in C. The drawback of this approach is 
that a complete triangulation of C is often too hard to accomplish. 

Instead of computing a full triangulation, we compute only a (regular) subdivision of C 
into few cones. To this end we remove one of the generators of the cones, say r s , compute 
the convex hull of the cone C = cone(ri, . . . , r s _i), and find all facets J 7 of C that are 
visible from r s . By J 7 ' we denote the set of all cones that we get as the convex hull of 
a facet in J 7 with the ray generated by r s . Then J 7 ' U {C} gives a regular subdivision 
of C, called the subdivision with distinguished generator r s . Before we now subdivide 
those cones in J 7 ' further into smaller cones, we use the following simple observation 
to remove cones that can be avoided due to the underlying symmetry given by S. 

Lemma 6. Let C,Cx, . . . , C/. C R™ be rational polyhedral cones such that C = Uf =1 Cj 
(not necessarily a disjoint union). Suppose that there is a permutation a and indices i 
and j such that Ci C o~(Cj) C C. Then the Hilbert basis of C is contained in the union 
of the Hilbert bases of the cones C\, . . . , Cj_i, cr(Cj), Ci+i, . . . , Ck- 

Proof. The result follows by observing that all lattice points in Cj also belong to o~(Cj) 
and thus can be written as a nonnegative integer linear combination of the Hilbert basis 
oiaiCj). □ 

If successful, this test whether Ci can be dropped is a very efficient way of removing 
unnecessary cones. However, the fewer generators are present in the cones Ci, . . . , Ck, 
the higher the chance that this test fails. So one has to make a trade-off between a 
simple test (that may fail more and more often) and a direct treatment of each cone 
Ci. As we compute only regular subdivisions whose cones are spanned by some of the 
vectors ri, . . . , r s , each of the cones C\, . . . , Ck can be represented by a characteristic 
0-1-vector x(C\), . . . , x(C&) of length s that encodes which of the generators of C are 
present in this cone. This makes the test Ci C o~(Cj) comparably cheap, as we only 
need to check whether x(Cj) < o~(x(Cj)). 

Summarizing these ideas, the symmetry exploiting approach can be stated as follows: 

(1) Let C = cone(ri, . . . , r s ) C R n and C = {C}. 

(2) i := 

(3) While C ^ do 

(a) i :— i + 1 

(b) For all K G C compute a subdivision with distinguished ith generator (if 
existent in K). 

(c) Let T be the set of all cones in these subdivisions. 

(d) Let M. be the set of those cones with a maximum number of rays. 

(e) Let C 7^ be the set M. together with all cones T G T that are not covered 
by a cone o~(M) with M G Ai and a G S, see Lemma El 
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(f) Remove from C all simplicial cones and compute their Hilbert bases. 

(4) For each computed Hilbert basis element h compute its full orbit {cr(h) : a G S} 
and collect them in a set %. 

(5) Remove the reducible elements from %. 

(6) Return the set of irreducible elements as the minimal Hilbert basis of C. 

This quite simple approach via triangulations and elimination of cones by symmetric 
covering already solves all three presented examples. In particular, it gives a compu- 
tational proof to Lemma [TJ The candidates for the representatives of Hilbert basis 
elements can be computed using "LattE for tea, too" by calling 

dest/bin/hilbert-f rom-rays-symm — hilbert-f rom-r ay s="dest/b in/hi lbert-f rom-rays" 

— dimension=26 S5.rays 
dest/bin/hilbert-f rom-rays-symm — hilbert-f rom-rays=" dest/bin/hilbert-f rom-rays" 

— dimension=43 355 . short .rays 
dest/bin/hilbert-f rom-rays-symm — hilbert-f rom-rays=" dest/bin/hilbert-f rom-rays" 

— dimension=42 346 . short . rays 



The data files can be found on http: //www. Iatte-4ti2 . de (For typographical reasons 
each command has been printed in two lines.) 

2.2. Second approach: partial triangulation. In the second approach, we build 
up a triangulation of the given cone C = cone(ri, . . . , r s ) C M n . However, by using the 
following Lemma [7] and its Corollary [HJ we can avoid triangulating many regions of the 
cone, since the triangulation would consist only of unimodular cones (for which the 
extreme ray generators already constitute a Hilbert basis), or, more precisely, avoid to 
construct simplicial cones whose non-extreme Hilbert basis elements are contained in 
previously computed simplicial cones. 

In the following we describe the facets of a full-dimensional rational cone by (uniquely 
determined) primitive integral exterior normal vectors. In other words, F = {x G C : 
c J x = 0} where c has coprime integer entries and c J y < for all y G C. 

Lemma 7. Let C = cone(ri, . . . , r&) C MJ 1 be a rational polyhedral cone such that 

• ri, . . . ,r k G Z n , 

• ri, ... , Tk-i Us in & facet of C defined by the hyperplane c T x = ; 

• c J r k = 1. 

Then the Hilbert basis of C is the union of {r^} and the Hilbert basis of cone(r 1; . . . , 
r fe -i). 

Proof. Let z G C n Z n . Then z = J2t=i X i r i for 

some nonnegative real numbers 

Ai, . . . , Afc. Multiplying by c T , we obtain 

k 

c T z = 22 X i° Jr i = AfcC T r fc = A fc . 
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As c, z G Z n , we obtain G Z. Hence, z is the sum of an nonnegative integer multiple 
of Tk and a lattice point z — A^r^. G cone(r 1; . . . , r^_i), which can be written as a 
nonnegative integer linear combination of elements from the Hilbert basis of this cone. 
The result now follows. □ 

This lemma implies the following fact, which excludes many regions when searching 
for missing Hilbert basis elements. 

Corollary 8. Let r%, . . . , G Z n such that C = cone(ri, . . . , r^-i) has dimension n, 
and C = C + cone(rfc). Suppose that r^ ^ C . Moreover, let F\,...,F q be the facets 
of C visible from and let ci, . . . , c q the normal vectors of these facets as introduced 
above. Then 

U{C) U {r fc } U |J {H(Fi + cone(r fe )) : |cjr fc | > 2, i = 1, . . . , g} 

generates C fl Z n . 

Proof. Evidently we obtain a system of generators of C fl Z n if we extend the union 
in the proposition over all facets F{, i — 1, . . . , q. It remains to observe that 

U(Fi + cone(r fc )) = {r k } U H{C D F<) 

if |cjrfc | = 1. But this is the statement of Lemma [71 □ 

Corollary [8] yields an extremely efficient computation of Hilbert bases — provided the 
case |cjrjfc| > 2 occurs only rarely, or, in other words, the system rj., . . . , of generators 
is not too far from a Hilbert basis. 

A thoroughly consequent application of Corollary [8] could be realized as follows, col- 
lecting the list A(C) of critical simplicial cones in a recursive algorithm. 

(CI) Initially A(C) is empty. 

(C2) One searches the lexicographically first linearly independent subset {r^,..., 
r,j d }. If the cone generated by these elements is not unimodular, it is added to 
A(C). 

(C3) Now the remaining elements among r 1; ...,r s (if any) are inserted into the 
algorithm in ascending order. Suppose that C is the cone generated by the 
elements processed already, and let r^ be the next element to be inserted. Then 
for all facets Fj of C such that cjr/j > 2 the list A(C) is augmented by A(Fi + 
cone(r i )). 

After all the critical simplicial cones have been collected, it remains to compute their 
Hilbert bases and to reduce their union globally, together with {ri, . . . , r s }. 

Let us add some remarks on this approach. 

(a) It is not hard to see that the list A(C) constitutes a subcomplex of the lexi- 
cographical triangulation obtained by inserting ri, . . . , r a . However, this fact is 
irrelevant for the computation of Hilbert bases. 
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4,392 


1,300 


# supp hyp 


4,948 


29,387 
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117,978 


# full tri 


2,654,000 


102,538,980 


? 


? 
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# partial tri 


48 


4,320 


775,800 


206,064 


3,109,495 


# cand 


96 


1,260 


41,593 


10,872 


168,014 



Table 1. Data of challenging Hilbert basis computations 



(b) In an optimal list of simplicial cones each candidate for the Hilbert basis of 
C would appear exactly once. (The candidates are the elements of the Hilbert 
bases of the simplicial cones.) The algorithm above cannot achieve this goal 
since the cones F + cone(r J ) are treated independently of each other. Neverthe- 
less it yields a reasonable approximation. 

(c) The drawback of the algorithm above is that it uses the Fourier-Motzkin elimi- 
nation recursively for subcones. Therefore Normaliz applies the algorithm above 
only on the top level and produces a full triangulation of the cones Fi + cone(rfc) 
for which c^r^ > 2 (instead of the list A(Fi + cone(rj))). 

(d) It is a crucial feature of the height 1 strategy that it reduces memory usage 
drastically. 

We illustrate the size of the computation and the gain of the improved algorithm by 
the data in Table [TJ In the table we use the following abbreviations: emb-dim is the 
dimension of the space in which the cone (or monoid) is embedded, dim denotes its 
dimension, # rays is the number of extreme rays, # HB is the number of elements 
in the Hilbert basis, # full tri is the number of simplicial cones in a full triangulation 
computed by Normaliz, # partial tri is the number of cones in the partial triangulation, 
# cand is the number of candidates for the Hilbert basis, and # supp hyp is the number 
of support hyperplanes. 

In addition to the improved algorithm just presented, parallelization has contributed 
substantially to the rather short computation times that (the experimental version of) 
Normaliz needs for the cones considered. The computations were done on a SUN Fire 
X4450 with 24 Xeon cores, but even on a single processor machine computation times 
would be moderate. 

Remark 9. Sturmfels and Sullivant [Tlj 3.7] stated a very interesting conjecture on 
the normality of cut monoids of graphs without i^-minors. For graphs with 7 and 8 
vertices we have used the approach via partial triangulations (and parallelization) in 
order to verify the conjecture. For these graphs no counterexample could be found. 
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